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Abstract 

We present a simple numerical method for the discrete breather construction based on the idea 
of the pair synchronization of the particles involved in the breather vibration. It can be used 
for obtaining exact breather solutions in nonlinear Hamiltonian lattices of different types. We 
illustrate the above method using chains of the coupled Duffing oscillators. With some additional 
approximation, the pair synchronization method leads to a very simple physical interpretation 
of the existence of the exact breathers as strictly time-periodic and spatially localized dynamical 
objects. 
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I. INTRODUCTION 



The concept of the discrete breathers or intrinsic localized modes in nonlinear Hamiltonian 

ri 

lattices have been proposed by Sievers and Takeno two decades earlier in Since that 
time, there appeared a vast number of papers considering different physical, mathematical 
and numerical aspects of the discrete breathers dynamics. The appropriate references can 
be found in many review papers (see, for example, {2], O, 0, [5]). 

Discrete breathers represent spatially localized and time-periodic excitations in nonlinear 
lattices. As we have already discussed in {f| , the possibility of the existence of such dynam- 
ical objects is not obvious. Indeed, because of spatial localization, different particles of the 
lattices vibrate with essentially different amplitudes. On the other hand, it is typical for 
nonlinear systems that frequencies depend on the amplitudes of vibrating particles. There- 
fore, one may ask: "How can a discrete breather exist as a strictly time-periodic dynamical 
object despite of its spatial localization" ? In the present paper, we will try to give an explicit 
answer to this question. 

There exist several rigorous existence proofs for breathers in networks of weakly coupled 
anharmonic oscillators reviewed, for example, in lal. The first existence proof of the dis- 

ri 

crete breathers was obtained by the principle of anticontinuity in [7j. The authors have 
demonstrated this approach for the chain of linear coupled nonlinear oscillators where anti- 
continual limit corresponds to zero coupling. The main idea of this approach is connected 
with application of the implicit function theorem in the space of time-periodic solutions at 
the discrete breather period. 

There exist a number of methods for breathers construction at computer accuracy pro- 
vided that their frequency (and its multiples) is far enough from resonance with linear 
phonon spectrum 0, 5]. As was emphasized by Aubry in 3], the Newton method is used 
practically in all such methods and "it is nothing but a numerical implementation of the 
implicit function theorem" used for rigorous proof of the DBs existence in Q. As a certain 
variation, in some numerical procedures, the Newton method is used for finding the Fourier 
coefficients of the exact breathers solutions 0, Q] . 

Unfortunately, the above methods, as well as the rigorous existence proofs, do not lead 
at least, from our point of view) to any simple physical interpretation of the DBs existence 

3)- 
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In the present paper, we develop a new simple method for breather construction (Sec. [TT]) 
which is able to produce breathers solutions with high precision (Sec. IHIj) and, on the other 
hand, leads to an explicit physical interpretation of DBs existence (Sec. HV|) . 

II. METHOD OF PAIR SYNCHRONIZATION 

Below we present the method of pair synchronization (PS-method) for discrete breather 
construction in nonlinear Hamiltonian lattices. Despite of the applicability of this method for 
wide class of such lattices, we prefer to demonstrate it with a simple case of the monoatomic 
chain of the linear coupled Duffing oscillators for clarity of exposition [15]. The generalization 
of the PS-method for other lattices seems to be more or less straightforward. 

The term "method of pair synchronization" is connected with the specific nature of its al- 
gorithm. Indeed, the PS-method represents an iterative procedure which implies subsequent 
synchronization of the vibrations of each pair of the breather's particles. 

The dynamical equations for the chain of the hard Duffing oscillators with linear coupling 
can be written as follows 

Xi + Xi + = j{x i+1 - 2xi + i = l..N. (1) 

Introducing the periodic boundary conditions 

x N+ i(t) = xt(t), x (t) = x N (t), (2) 

we will consider the simple symmetric breathers in the chains (1) with iV = 3, 5, 7, 9. ..etc. 
particles. Because of the periodic boundary conditions (121), we can imagine that our chain 
represents a ring depicted in its equilibrium state in Figs. [Q 0, [T^J 

We start with the case N = 3 (see Fig. [I]). Let the first particle be the central particle 
of the breather and, therefore, its amplitude is greater than amplitudes of the peripheral 
particles, i.e. 

\xi(t)\>\x 2 (t)\,\x 3 (t)\. (3) 

Moreover, the relation x 2 (t) = x 3{t) must hold as far as we search for the symmetric breather. 

To construct such breather, the following initial conditions can be used for integrating 
the dynamical equations ([!]): 

zi(0) = x 2 (0) = x 3 (0) = 0, Xl {0) = at, x 2 (0) = x 3 (0) = a 2 . (4) 
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FIG. 1: Scheme of the initial conditions for symmetric breather in the chain ([T]) with N = 3 
particles and periodic boundary conditions: a?i(0) = X2(0) = £3(0) = 0; xi(0) = a±, £2(0) = 02, 
±3(0) = a 3 . 




FIG. 2: Start of the pair synchronization method for the chain with N = 3 particles [xi(0) = 2, 
i 2 (0) =x 3 (0) =0]. 




FIG. 3: The functions xi(t) and X2(t) for different values of i;2(0) = «2 for the model (pQ) with 
7 = 0.3 [4 0) = -0.7, 4 1} = -°- 4 ) 4 2) = 

Thus, we suppose that all particles of the chain are in their equilibrium positions at the 
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FIG. 4: Duffing chain ([!]) with N = 3 and 7 = 0.1: the final result of the synchronization of the 
vibrations of the particles {et\ = 2, a 2 = —0.144656). 

initial instant, but they possess certain velocities which we denote by a i: (i = 1..JV). 

There are two parameters from which our breathers depend on, namely, 7 entering the dy- 
namical equations ([T]) and the initial velocity a,\ of the central breather's particle. Choosing 
some reasonable values of these parameters and supposing a 2 = 0, we integrate equations (pQ) 
for N = 3 over time interval containing, at least, one zero of the function X\{t) and one zero 
of x 2 (t) (see Fig. [2]). Now let us vary the initial velocity x 2 (0) = a 2 of the second particles 
aiming at coincidence of the first zero of x 2 (t) with that of x\{t). In other words, we try to 
synchronize the oscillations of the first particle with those of two neighboring particles 

[x 2 (t) =x 3 (t)]. 

Different numerical methods can be applied for such purpose, for example, the shooting 
method using the idea of dichotomy. Zeros of the functions Xi(t) and x 2 (t) can be found by 
the Newton-Raphson method or by some other well-known numerical methods. 

Some comments are appropriate at this stage. 

Let i\ = 1, 2, 3, ... labels the subsequent zeros of the function Xi(t) and i 2 = 1, 2, 3, ... 
labels those of the function x 2 (t). In principle, we can superpose zeros of the functions xi(t) 
and x 2 {t) with different values i\ and i 2 . As a consequence, we obtain discrete breathers of 
different types (see Sec. IHIj) . 

In this section, we consider only the single-frequency breathers, i.e. the breathers for 
which all particles of the chain oscillate with one and the same frequency ujj = ooj, (j = 1..JV). 
This is not a general case. Indeed, one can imagine two-frequency, three-frequency, etc. 
breathers whose particles possess different, but multiple frequencies (for example, if u)\ = 2u 2 
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or uj\ = 3u>2 for the chain with N = 3 particles, we have two- frequency breather). In such 
a case, the breather frequency Uf, is equal to the minimal frequency among all possible 
Uj (j = 1..N) of the individual particles. Therefore, all the particles vibrate with the 
frequencies which prove to be multiples of cj&. 

In order to avoid some misunderstanding, let us note that the single-frequency breather 
is not monochromatic, i.e. although all the breather's particles possess identical vibrational 
frequencies, the Fourier series of periodic functions Xi(t) contains not one, but a certain set 
of different harmonics. 

Some comments on the terminology are appropriate at this point. The terms single- 
frequency and many-frequency breathers, which we use in the present paper, disagree with 
those used in jg]. Indeed, in the latter paper, a single-frequency discrete breather means 
the conventional strictly time-periodic dynamical object, while a many-frequency breather 
discribes a quasiperiodic motion on an n-dimensional torus (n basic frequencies correspond 
to its Fourier decomposition). In contrast, considering only strictly time periodic breathers, 
we divide them into different classes: in single-frequency breathers, all the particles vibrate 
with one and the same frequency, while in many-frequency breathers, different particles 
possess different, but divisible frequencies (see below for more details). 

Let us return to the discussion of the numerical procedure of the pair synchronization for 
the chain with N = 3 particles. 

In Fig. [3] we depict the function x^ify for several values a.2 of the initial velocity of the 
second particle. From this figure, one can see how the position of the first zero of X2(t) 
changes with changing £2(0) = 0.2- Note that changing of 02 leads not only to a shift of the 
zero of X2(t), but also to a certain shift of Xi(t), but the latter turns out to be rather small 
and is not noticeable in Fig. [31 

It is obvious from Fig. [3j that oup and provide a "fork" to which we can apply the 
method of dichotomy to attain concidence of the zeros of the functions X\{t) and Xzif) with 
a chosen level of accuracy. The final result of the discussed synchronization is shown in 
Fig. HI where one can see a single-frequency breather for the chain with N = 3 particles. 

The synchronization of the vibrations of the first and second particles we call 5(1, 2) pro- 
cedure. Complete realization of this procedure for the case of the single-frequency breather 
leads to the full coincidence of all the zeros of the functions Xi(t) and ^(i) with equal labels 
%\ and %i. In general case, the similar procedure of synchronization of the vibrations of i-th 
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FIG. 5: Scheme of the initial conditions for the five-particles chain. 

and j-th particles will be called S(i,j) procedure. 

Now, let us construct a single-frequency symmetric breather for the chain (CQ) with N = 5 
particles (see Fig. [5]). Choosing resonable values of 7 and aci, we have to find such values of 
ct2 and a?3 which lead to the equal vibrational frequencies of all the particles of our chain. 
In other words, we want to synchronize vibrations of the central particle with those 

of peripheral particles [x 2 (t) = x 5 (t) and x 3 (t) = x±(t)]. 

Let us apply the method of pair synchronization for discrete breather construction in the 
considered chain. The algorithm of the method can be described as follows. We syn- 
chronize vibrations of the particles 1 and 2 with aid of the above described procedure 

n 

5(l,2)[17j]. Then we synchronize vibrations of the particles 2 and 3 with the aid of the 
procedure 5(2,3). The sequence of these two procedures, 5(1,2), 5(2,3), forms one it- 
erative cycle of the pair synchronization method. We must repeat such cycles, namely, 
5(1,2); 5(2,3)| 5(1,2); 5(2, 3)| 5(1,2); 5(2, 3)|..., up to attain a desirable level of ac- 
curacy. Indeed, the coincidence of zeros of xi(t) and X2(t), after finishing the procedure 
5(1,2), will be slightly disturbed because of applying the procedure 5(2,3) which changes 
0:2, CK3, and we must recorrect the accuracy of this coincidence by the next iterative cycle. 
One can hope for convergence of such iterative process, at least, for the case of sufficiently 
strong space localization of the discrete breather. 

The generalization of the considered method to a chain with arbitrary iV is straightfor- 
ward: we repeat iterative cycles 

{5(1,2); 5(2,3); 5(3, 4); ... S(N - 1, N)} (5) 

up to obtaining the discrete breather with resonable precision. Note that the method of 
pair synchronization can be regarded as a version of the relaxation procedure used in many 
different numerical methods. 
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Note that considering chains with subsequently increasing number of particles (N = 
3, 5, 7, 9, 11, ...), we can use the final set of initial velocities (ai, «2, «3,... ) of the chain 
with iV = m in order to begin calculations for the chain with N = m + 2, assuming two last 
velocities equal to zero. 

III. SOME NUMERICAL RESULTS 

Below, we present some results of the discrete breathers construction with the aid of the 
pair synchronization method for a number of mechanical models and for different kinds of 
discrete breathers. The application of this method to other oscillatory chains and multidi- 
mensional lattices will be discussed in future publications. 

A. Linear coupled Duffing oscillatory chain 

Let us consider the model ([1]) for different values of the coupling parameter 7 and different 
values of the initial velocity ati of the central breather's particle. 

The appropriate initial velocities x«(0) = <Xj, i = 2..N are necessary for a start of the pair 
synchronization procedure. For 7 1 one can obtain these velocities from the approximate 
formula, which is derived in Sec. HVl 

a i+1 = 2 2 a " (i = l-N-l). (6) 

u b - u 

Here Ub is the breather frequency, which can be roughly approximated as Ub = ou(ati) with 
ou(oii) being the frequency of a single Duffing oscillator, while Uq = 1 + 27. The function 
u>(a) is determined by the solution of the Duffing equation x + x + x 3 = and we depict it 
in Fig. [6] for obtaining rough estimates of the breather frequency uo^. 

In Table [H we present the velocity and amplitude breather profiles, as well as the breather 
frequency Ub, for the chains with different number of particles iV = 3, 5, 7, 9, 11, 13 for 
7 = 0.3, cui = 2. In the first column of each fragment of this table the indices (n) of 
the particles are given (for the central breather particle n = 0). In the following columns, 
we present the initial velocity (a n ) and vibrational amplitude (A n ) of each particle of the 
considered chain. 

From this table, one can see how the discrete breather for the infinite chain acquires its 
shape step-by-step with increasing N. We restrict our consideration to N = 13 because, 
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FIG. 6: The function to (a) determining dependence of the Duffing oscillator frequency uo on the 
initial velocity a. 




FIG. 7: The amplitudes A n for the symmetric breathers in the Duffing chain (pQ) for different values 
7 (in the logarithmic scale). 

in the case 7 = 0.3, the amplitudes of the particles which are the most distant from the 
breather's center do not exceed 10~ 4 of the amplitude of the central particle. 

In Table [TTJ we present information similar to that in Table [I] for the chain with N = 9 
particles, but for different values of the coupling parameter 7 (0.3 < 7 < 1.5). From this 
table, one can see that a degree of localization becomes worser with increasing 7. 

Because the amplitudes A n of the particles decrease rapidly with increasing their distance 
from the breather center, we depict them in the logarithmic scale for different 7 in Fig. [7J 
From this picture one can see the exponential decay of the amplitudes of the peripheral 
particles. The degree of the breather space localization decreases with increasing 7. 

Let us note that the pair synchronization method can be often used not only for suffi- 
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TABLE I: Discrete breathers for the Duffing chain ([!]) for 7 = 0.3 and different number of particles 
(AO. 
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ciently small coupling between the oscillators of the lattice. From Table [TTJ we see that 
this method also works, in the case of the model ([I]), for 7 ~ 1. Moreover, the pair 
synchronization method can be used for 7 > 1 and even for 7 ^> 1 (in such cases for- 
mula ([6]) cannot be used for obtaining an initial velocity profile). As an example, let 
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TABLE II: Discrete breathers for the Duffing chain (TT]) for different values of the coupling parameter 
7 (N = 9). 
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us give the set of the initial velocities x«(0) = ctj obtained by this method for the 
single-frequency discrete breather in the chain ([1]) with N = 7 particles for 7 = 20: 
{«i = 2, a 2 = -1.801427, a 3 = 1.245828, a 4 = -0.444403}. Existence and stability 
of the discrete breathers in chains with large coupling will be discussed elsewhere. 
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FIG. 8: Two-frequency symmetric breather in the Duffing chain ([!]) for 7 = 0.3, ol\ = 7, N = 5. 
B. Many-frequency discrete breathers 

First of all let us comment on one point of the procedure of pair synchronization. Consid- 
ering this procedure in Sec. HIl we have aimed at coincidence of those zeros of two functions, 
Xi(t) and Xj(t), which are most close to the initial instant t = 0. This is not a general 
case. Indeed, one can imagine that Xi(t) and Xj(t) possess equal periods T{ = Tj = T, but 
that they do not go to zero simultaneously at any internal instant < t < T. Then we 
should aim at coincidence of the second zeros of the considered functions, i.e. the relation 
Xi(T) = Xj(T) = must hold. 

Now let us focus on constructing many-frequency discrete breathers. Let us remember 
that we coin this term for such spatially localized dynamical objects whose particles vibrate 
with several different but divisible frequencies. This means that many-frequency breather 
proves to be a time-periodic object whose period T is the largest of the multiple periods 7\, 
T 2 , T/v of all the particles of the chain. 

In Fig. [8], we depict a two-frequency symmetric breather (00% = 2uj2 = 2003) in the chain 
(CQ) with N = 5 particles for 7 = 0.3. Indeed, the first particle vibrates with the 

frequency twice larger than those of all other particles [^(i), x 3(t)]- This breather with 

= ^2(= ^3) = 1-39615 corresponds to the following data: 7 = 0.3, 07 = 7, a 2 = 
— 1.064224, as = 0.352247. Hereafter, we do not point out explicitly that all the initial 
displacements are zero: Xj(0) = 0, i = 1..N. 

Note that for the same 7 and 07, we can find a single-frequency breather (oj-y = U2 = 
U3 = 2.77828), as well (see Fig. [9] which was obtained for 7 = 0.3, a% = 7, a.2 = —0.379206, 
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FIG. 9: Single-frequency symmetric breather in the Duffing chain (TTJ) for 7 = 0.3, ol\ = 7, N = 5. 
Vibrations of the third particle is not visible, because of its small amplitude. 




FIG. 10: Two-frequency symmetric breather in the Duffing chain ([I]) for 7 = 0.3, a± = 19, N = 5 
(uji = 3u>2 = 3^3). This breather is unstable. 

a 3 = 0.017876, N = 5). 

In Fig. dni we depict the two-frequency breather with Uf, — U2 — 1.48666 whose particles 
possess the following frequencies: uj\ = 3u; 2 = 3u 3 . This breather corresponds to the 
following initial data: 7 = 0.3, at = 19, a 2 = -1.593445, a 3 = 0.456080, N = 5. 

In Fig. [TlJ a three-frequency breather with U\ = 2c<j 2 , c<j 2 = 2u>3 is shown. It corresponds 
to the following initial data: 7 = 0.3, on = 28, a 2 = -7.018426, a 3 = 1.472607, = 5. 
The breather frequency is uji, = u 3 = 1.35108. 

Note that the pair synchronization method providing us with an initial breather profile 
does not guarantee the breather stability. For example, the breather depicted in Figs. dUJ dU 
proves to be unstable dynamical object, and this fact can be checked by straightforward 
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K(t) 




FIG. 11: Three-frequency symmetric breather in the Duffing chain ([I]) for 7 = 0.3, a\ = 28, 
N = 5 {uji = 2u)2 = 4u;3). This breather is unstable. 




FIG. 12: A nontrivial two-frequency breather in the Duffing chain (pQ) for 7 = 0.3, a\ = 4, N = 5 
[iOl = 2u>2 = 2u>3 = 2to^). 

integration of the dynamical equation ([T|), as well as with the aid of the Floquet method. 

Finally, in Fig. [T2J, we demonstrate a nontrivial example of two-frequency breathers (7 = 
0.3, a ± = 4, a 2 = -0.233618, a 3 = 0.350414, N = 5). Indeed the function x 2 (t), unlike x 3 (t), 
possesses two different maxima (minima) within one breather period which coincides with 
that of X2(t) and ^(t). Note that vibrational amplitudes of the second and third particles 
are of the same order of smallness. 

Note that, for simplicity, in many examples of this section we consider only chains with 
small number of particles (iV). This is not an essential restriction because, in every such 
case, the synchronization procedure can be easily continued for much greater values of N. 

Let us emphasize once more that, unlike the terminology of the paper [8], our many- 
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FIG. 13: Single-frequency antisymmetric breather in the Duffing chain ([!]) for N = 6, 7 = 0.3, 
ai = 2. 

frequency breathers are time-periodic dynamical objects, but different particles vibrate with 
different (divisible) frequencies Ui (i = 1..N). We would like to note that one can found in ^ 
a certain mention of multibreathers whose particles vibrate with different but commensurate 
frequencies. However, we do not know any papers with detailed analysis of such dynamical 
objects and it seems that our many-frequency breathers are not fully identical with the 
above mentioned multibreathers. 



C. Constructing of antisymmetric breathers 



Up to this point, we have considered only the symmetric discrete breathers, i.e. breathers 
of the form {...x 3 (£), x^it), x\(t), xo(t), Xi(t), X2(t), xs(t), ...}. Often they are called Sievers- 
Takeno modes. 

The pair synchronization method can be also used for constructing antisymmetric 
breathers which possess the following form {... — x$(t), —X2(t), — xi(t), x±(t), X2(t), x^(t), ...}. 
Such breathers are usually called Page modes. In Table llTlj we present initial velocities (a*) 
and amplitudes (Ai) of all particles of the antisymmetric breather for the Duffing chain (pQ) 
with N = 8 particles. One can see the time evolution of this breather in Fig. [13j 
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TABLE III: A single-frequency antisymmetric breather in the Duffing chain ([I]) for different values 
of the coupling parameter 7 (N = 8). 
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D. Other types of oscillatory chains 



Chains of the coupled Duffing oscillators are used for some physical applications. For 
example in 9|, the model of such a type is exploited for describing discrete breathers in 
cantilever arrays. Indeed, the authors of that paper consider an oscillatory chain with on-site 
and inter-site forces both containing linear and cubic terms. Now we want to demonstrate 
that the pair synchronization method can be applied for the breather construction in the 
different chains of such class of mechanical systems. 
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We will consider the following two models: 

it + Xi + x\ = 7[(x i+ i - xif - (xi - ^-i) 3 ], (7) 

and 

Xi + x\ = j[(x i+1 - XiY - (xi - Xi_i) 3 ], (8) 

with periodic boundary conditions (j5]). 

The model (J7J) contrary to (pQ) represents the chain of the Duffing oscillators with cubic 
coupling. The model (|SJ) represents the so called chain, i.e. the chain whose potential 
energy is a uniform function of the fourth order. Discrete breathers in the K4 chains with 
and without on-site interactions were studied in many papers from different points of view 



0,Q, 



111 ]. This chain is very convenient for investigation because it allows separation of 
space and time variables. In our paper [fjj], the symmetry discrete breathers were found for 
the model (jHJ) with high precision using the concept of the similar nonlinear normal modes 



presented in jfj]. 

In Table IIV| we present the velocity and amplitude profiles for the chain (J7|) for some 
values of the coupling parameter 7 in the chain with N = 7 particles. It is worth comparing 
this table with that for the linear coupled Duffing oscillators (see Table ITT]) . One can see that 
the replacing of the linear coupling of the model (CQ) on the cubic coupling of the model (|7j) 
changes the degree of the breather spatial localization. 

In Table |Vj we present the similar information for the K4 chain, i.e. for the chain with 
cubic force interactions and without the linear on-site forces. As was already mentioned, 
for the latter model there exist Rosenberg nonlinear normal modes. For these modes, the 
displacement of every particle of the chain at any instant t is proportional to the displacement 
of an arbitrary chosen particle, say, the first particle. Therefore, for the dynamical regime 
described by a Rosenberg mode one can write 

x i (t) = k i -x 1 (t), i = l..N, (9) 

where hi are the constant coefficients which can be determined from a certain systems of 
nonlinear algebraic equations (see 6j for the further details). 

From ([9]), one can see that Xj{t) = fc i x 1 (t), i = 1..N and, therefore, Xj(0) = fcjii(O). In 
turn, this means that the initial velocity profile is proportional to the spatial profile of the 
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introduced by Rosenberg in [12|. The detailed investigation of the breather stability is also 



TABLE IV: Symmetric discrete breathers for the chain (JT]) for different values of the coupling 
parameter 7 (N = 7). 
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Rosenberg mode which was found for some range of 7 in [fij]. The results of that paper and 
those from Table |V] of the present paper, obtained by the method of pair synchronization, 
proved to be identical. 

IV. A SIMPLE PHYSICAL INTERPRETATION OF THE PAIR SYNCHRONIZA- 
TION PROCEDURE 

In Sec. Ull we have described the method of pair synchronization for discrete breather con- 
struction. At each step of this method, we force two particles, who vibrate with essentially 
different amplitudes, to have exactly identical frequencies. Keeping in mind that in nonlinear 
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TABLE V: Symmetric discrete breathers for the chain (|8|) for different values of the coupling 
parameter 7 (N = 5). 



1 


= 0.3; cob - 


= 1.86714 


0'' 
/ 


= 0.5; = 


= 1.97301 


l 


= 0.7: tOh = 


= 2.13965 


n 


iJL n 


A 


11 


rv 
u n 


A 


n 




A 


-2 


0.0071987 


0043733 


_2 


017014 


010332 


_2 


023577 


013202 


-1 


-0.598577 


363639 


-1 


-0 776917 


471793 


-1 


-0 857977 


0.480441 





2 


1.215014 





2 


1.214528 





2 


1.119939 


1 


-0.598577 


0.363639 


1 


-0.776917 


0.471793 


1 


-0.857977 


0.480441 


2 


0.0071987 


0.0043733 


2 


0.017014 


0.010332 


2 


0.023577 


0.013202 


7 = 1; u b = 


2.33466 


1 


= 1.2; uj h = 


= 2.44172 


'Y 


= 1.5: ujh - 


= 2.57983 


n 


«n 


An 


n 






n 




An 


-2 


0.029452 


0.015115 


-2 


0.031948 


0.981391 


-2 


0.034559 


0.016050 


-1 


-0.917511 


0.470863 


-1 


-0.940078 


0.461292 


-1 


-0.962237 


0.446887 





2 


1.026392 





2 


0.981391 





2 


0.928850 


1 


-0.917511 


0.470863 


1 


-0.940078 


0.461292 


1 


-0.962237 


0.446887 


2 


0.029452 


0.015115 


2 


0.031948 


0.981391 


2 


0.034559 


0.016050 



regimes, particles vibrating with different amplitudes possess, as a rule, different frequencies, 
one can wonder what is the physical nature of the possibility of the above discussed pair 
synchronization? The understanding of this nature leads us to the explicit answer to the 
question formulated in the introduction of the present paper: "How can discrete breathers 
exist as space localized and, at the same time, strictly time-periodic dynamical objects"? 

Let us consider the chain ([1]) for N = 3 and write down explicitly the corresponding 
dynamical equations: 

xi + %i + x\ = 7(22 - 2a;i + x 3 ), 

x 2 + x 2 + x\ = 7(x 3 - 2x 2 + xx), (10) 

x 3 + x 3 + xl = 7(0:1 - 2x 3 + x 2 ). 

These equations describe the vibrations of three Duffing oscillators with linear coupling 
and parameter 7 determines the strength of this coupling. The relation x 2 (t) = x 3 (t) holds, 
since we search for a symmetric breather. Therefore, tree equations ( fTUj) reduce to the 
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following two equations: 

Xx + X\ +xl — 2j(x 2 — Xi), (11) 
x 2 + x 2 + xl = 7(xi - x 2 ). (12) 

We suppose 

0<7<1, \x 2 (t)\ < (13) 

Actually, the last relation is the condition of strong space localization of our three-particle 
breather. Then it is resonable to neglect the term x\ in Eq. (fT2|) and the term 27x2 in 
Eq. ffTTj) . As a consequence of such simplification, Eqs. (Illtil2p read [18j 



xi + (1 + 2 7 )x! + x\ = 0, (14) 

x 2 + (1 + t)x 2 = 7x1. (15) 

Eq. (|T4|) represents the Duffing equation whose solution can be expressed via Jacobi elliptic 
functions, while Eq. ffT5l) is the equation of the harmonic oscillator with periodic external 
force 7 • Xi(t). 

Let us make a next step of simplification replacing periodic function X\(t) in ffTBI) by its 
first Fourier harmonic or, more precisely, let us assume 

Xi(t) = Asin(ujt) + Bcos(ujt), (16) 

where the correct frequency u can be obtained as a result of solving the Duffing equation f|T4l) . 
From the initial conditions 

xi(0)=0, xi(0) = «i, (17) 
we find for the solution (jTBT) that A = — , 5 = and, therefore, 

xi(t) = — sin(cut). (18) 

In Fig. [TH we depict the exact solution to Eq. (fl4"|) (by a solid line) and the approximate 
solution ([IB"]) (by a dotted line) for 7 = 0.1, a\ — 2. From this figure, one can see how 
close can be to each other the above two solutions. 

Substituting f[TB"j) into r.h.s. of Eq. ffTBI) . we can write 

x 2 (t) + uj ■ x 2 = sm(cut), (19) 
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FIG. 14: The exact solution to the Duffing equation (|14p for 7 = 0.1, ol\ = 2 (solid line) and its 
approximation (fT8|) (dotted line). 



where c<Jq = 1 + 7. Eq. ( |T9i) is the equation of the harmonic oscillator with fundamental 
frequency ujq and with periodic driven force proportional to sin(a>£). 

Taking into account the initial conditions £2(0) = 0, £2(0) = «2, we easily obtain the 
solution to Eq. ffT9]) : 



£2(2) 



7«i 

«2 



^0 - 



sin(d;ot) H r^? ^sin(cut). (20) 

a;(o;Q — ur) 



Here the first term containing sin(d;o^) represents a solution to the equation ffl9l) without 
right-hand side, i.e. a solution of the free harmonic oscillator with frequency uq, while the 
second term containing sin(u;t) represents a particular solution to the full Eq. ffl9|) with the 
periodic driven force in its r.h.s. 

Let us analyze the solution (|2"U|) . In general case, the frequencies ujq and ui are incommen- 
surable. Therefore, for existence of the discrete breathers as strictly time-periodic dynamical 
objects it is necessary to eliminate the contribution from the "vibrational individuality" of 
the peripheral particles, i.e. the coefficient in front of sin(u;ot) must be equal to zero. From 
this condition, we find the initial velocity a<i of the second particle via that («i) of the 
central particle: 

a 2 = 5 2. (21) 

Note that uo > ujq since the frequency of the hard Duffing oscillator increases with increasing 
of the vibrational amplitude. Therefore, a.2 and cti possess opposite signs. 

Certainly, Eq. (I2T!) gives us only a certain approximation to «2, because we take into 
account only one harmonic from the Fourier decomposition of the function X\ (t) . Thus, our 
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approximate breather solution is 



xt(t) = — ■ sm(ut), (22) 



x 2 (t) = x 3 (t) = — 2T ■ sin (art). (23) 

LO\UJ A — UJq) 

It is essential, that the vibrational amplitude of X2(t) is ^/{oj 2 — Uq) times less than the 
amplitude of the central particle. In a certain sense, the relation 7 C 1 means the "space 
localization" of our breather in the three-particle chain. 

Actually, obtaining (I22H23P we have applied the synchronization procedure 5(1, 2) (see 
Sec. (TT]) in the framework of our approximation to nonlinear dynamics of the Duffing chain ([1]) 
with N = 3. Considering the chain with N = 5 particles (see Fig. [5]), we must also 
synchronize the vibrations of the particles 2 and 3. Obviously, all arguments which brought 
us to Eq. (12211231) can be applied to this synchronizing. Indeed, because of Eq. (j2"3~l) . the 
second particle x^if) = x^(t) already vibrates with the frequency u of the central particle. 
Therefore, the procedure 5(2, 3) reduces to the finding solution for x 3 (t) = x 4 (t) by solving 
the harmonic oscillator equation for x 3 (t) with time-periodic force driving exerted by the 
function x%{t) from Eq. ( f2~3l) . Demanding full suppression of the fundamental frequency 
contribution of the third particle [function x$(t)\ leads us to the relation similar to (12T1) 



«3 = (24) 



W — LU, 







where Uq = 1 + 7. 

Thus, the vibrational amplitudes of the third and fourth particles -jp^ji times less than 
that of the particles 2 and 5 (see Fig [5]). 

Obviously, all this argumentation can be continued to the consideration of the breathers 
in the chains with arbitrary number of particles. As a result, we conclude that vibrational 
amplitudes of all the peripheral particles decrease, at least, by the factor 7 -C 1. In turn, 
this means that our discrete breather possesses exponential decay in vibrational amplitudes 
when we pass on from the central particle to the more distant (peripheral) particles. 



One additional comment, concerning the fundamental frequencies ojq = y/1 + 27 and 
cJ = \/\ + 7 is appropriate at this point. 

For the symmetric breather in the chain with N = 5 particles (see Fig. [5]), there hold 
relations Xz{t) = xs(t), x 3 (t) = x^(t). Substituting these relations into the dynamical 
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FIG. 15: Scheme of the initial conditions for the seven-particle chain, 
equations (1), we obtain: 

X\ + Xi + x\ = 27(22 - Xl), 

x 2 + x 2 + xl = 7(23 - 2x 2 + xi), (25) 

x 3 + x 3 + x\ = 7(22 - x 3 ). 

Then taking into account the relations |&2(0I ^ l x i(^)l an d l x 3(^)l ^ I ^2^)1; which are a 
consequence of the breather localization, one can reduce Eqs. fl25l) to the following form: 

'±\ + uj^xi + xl = 0, 

x 2 + u\x 2 = 7^1, (26) 

x 3 +u>%x 3 = 7x2, 

containing slightly different fundamental frequencies ujq and ujq. 

The similar procedure for the case N = 7 with taking into consideration the relations 
(see Fig. [TSj) x 2 (t) = 27 (t), x 3 (t) = xe(t), x<i(t) = x$(t), brings us to the equations 

Xi + UJqXi + x\ = 0, 
X 2 + UqX 2 = IXx, ^ 
X S +U%X 3 = 7x2, 
X 4 + UqX 4 = 7x3. 

Proceeding in such a way for the chains with greater N, we conclude that all approximate 
equations contain fundamental frequency u>o = y/1 + 27 excepting only the last equation 
which contains the fundamental frequency uq = y/1 + 7. Bearing in mind that 7 < 1, we 
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obtain, from the equations similar to Eqs. (I26I) and (I2T1) . the above formulated conclusion 
about exponential decay of the amplitudes of the peripheral breather's particles. 

Let us summurize the above results. It has been manifested that the existence of discrete 
breather is connected with the "suppressing of the individuality" of all peripheral particles by 
vibrations of the central particle of our breather. In other words, the central particle forces 
upon the peripheral particles its own rhythm of vibrations. To obtain the correct discrete 
breather one must tune onto this rhythm with high precision by suppressing all terms in the 
solution to Eq. ([1]) corresponding to vibrations with different frequencies of all the peripheral 
particles. Therefore, an exact discrete breather seems to be a very unusual dynamical object 
and, in practice, we deal with quasibreathers rather than with exact breathers 6|. For 
the former objects, it is typical that all breathers particles vibrate with slightly different 
frequencies and these frequencies slightly drift in time. In turn, this means that in the 
spatially localized solution to Eq. ([1]) there exist terms with slightly different fundamental 
frequencies of the peripheral particles, whose amplitudes are sufficiently small. 

We believe that the presence of small amplitude terms with the frequencies ujj different 
from the principal breather frequency Ub, peculiar to any quasibreather solutions, cannot 
be the cause of the loss of their stability (see, for example, Eq. (1201) ). In other words, the 
loss of the strict periodicity does not mean, in general, the loss of the spatial localization of 
the considered dynamical objects. However, at the present time, we cannot give a rigorous 
mathematical proof of this proposition. 

In conclusion, let us make two additional comments. 

1. For the chains with an uniform potential, one can construct the Rosenberg nonlinear 
normal modes 6j. Our chain ([T|) of the Duffing oscillators does not belong to such 
class of mechanical systems. Nevertheless, the approximate solutions obtained in this 
section [see, for example, Eqs. (I22H23P ] prove to be the Rosenberg mode, since ratios 

= kj (j = 1..N) do not depend on t. But if we take into account the next 
Fourier harmonics for Xi(t), such ratios do depend on time. We illustrate this fact in 
Fig. HE] for the chain with N = 3 particles. 

2. In previous sections, we consider not only the single-frequency breathers, but many- 
frequencies breathers, as well. The existence of such dynamical objects can be also 
understood in the framework of our approximate methods for the discrete breathers 
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FIG. 16: The ratio for the exact solution to Eqs. ([TTj IT2 ]) for -y = 0.3. 



construction discussed in the present section. Indeed, let us return to the expres- 
sion (I20p for the chain with N = 3 particles. Now we will not annihilate the coefficient 
in front of sin(ujQt). Instead of this, we choose such value of the initial velocity of the 
central particle that leads to the frequency uj, which is a certain multiple with respect 
to the frequency ujq of the peripheral particle: uj = muo with m being an integer 



number 



19j. Then the resulting breather possesses the frequency ujq, determined by 



the vibrations of two peripheral particles (2 and 3), while the frequency u = mujQ of 
the central particle turns out to be higher multiple of ujq. Thus, we obtain a discrete 
breather whose particles vibrate with two different, but divisible, frequencies (ou and 
uj = mujo), i-e. we find a two- frequency breather. In this case, the strictly periodic 
discrete breather exists also only if we tune onto the relation uj = mujQ with high 
precision, otherwise we obtain a quasibreather with two incommensurable frequencies 
uj and ujq. 



When this paper was finished, we became aware of the very interesting paper [lc 
by Ovchinnikov concerning the localization of the vibrational energy in strongly excited 
molecules and molecular crystals. The author has presented the elegant arguments for the 
possibility of the above localization using as an example a molecular dimer described by 
equations similar to Eqs. (11 1111 21) of our paper. The everaging procedure was applied to 
obtain an approximate analytical solution to these dynamical equations in order to demon- 
strate the existence of the periodic solution with different vibrational amplitudes of the 
dimer particles. Above in this section, we have presented somewhat different arguments for 
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the same goal. 



V. CONCLUSION 

In this paper, we present a simple method for the discrete breather construction which has 
been called the pair synchronization method because of its transparent physical meaning. 
Indeed, it represents a certain iterative procedure which, at each step, synchronizes subse- 
quently the motion of the pairs of particles involving in the breather vibration. We believe, 
that this method can be applied for construction exact breathers in the nonlinear lattices 
of different types. The Duffing oscillatory chains with linear and nonlinear coupling were 
considered for explanation how the pair synchronization method can be used in practice. 

A comparison of efficiency of the pair synchronization method with other methods for the 
discrete breather construction based on the Newton procedure for solving nonlinear algebraic 
equations, on the steepest descent procedure etc. will be presented elsewhere. 

With some additional approximation, the pair synchronization method leads to a very 
simple physical interpretation of the possibility of the exact breathers existence as strictly 
time-periodic and spatially localized dynamical objects in nonlinear Hamiltonian lattices. 
Indeed, constructing a strictly time-periodic breather, we actually demand annihilation of 
the contributions with natural vibrational frequencies from all the peripheral particles to 
the breather solution. This means the existence of some kind of dictatorship of the central 
particle (in general, of the breather core) which must "suppress any individuality" of the 
peripheral particles by above mentioned annihilation. In other words, all terms of the exact 
breather solution with frequencies different from that of the central particle (which appear 
from the natural frequencies of the peripheral particles) must be turn into zero in compulsory 
order). The breather core compels the peripheral particles vibrate with its own frequency. 
In Sec. IIV} we demonstrate, for the case of weakly coupled oscillatory chains, that this 
compulsion leads to the exponential decay of the vibrational amplitudes of the peripheral 
breather particles. 

On the other hand, if we do not tune the initial conditions on the exact breather solution, 
i.e. if the vibrational contributions from the peripheral particles are not equal to zero, we 



obtain a certain quasibreather [6]. This spatially localized dynamical object is characterized 
by different frequencies appearing from the peripheral particles. The difference in frequen- 
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FIG. 17: An extremely long-lived quasibreather for the Duffing chain ([I]) with N = 5 particles 
which, possibly, is a stable dynamical object [7 = 0.3, a\ = 2, 02 = —0.5760, 03 = 0.1]. 



cies is brought about by the phenomenon typical for nonlinear dynamics, namely, by the 
dependence of the frequency on the vibrational amplitude. 

In (fl, we proposed to characterize the proximity of the quasibreather to an exact time- 
periodic breather with the aid of the mean square deviations in frequencies of the individual 
particles and that of the fixed particle in time (there exist a certain temporal drift of the 
frequency of every chosen particle). We also presented there some arguments (at least, for 
the case of K4 chain) for the possibility of the quasibreather stability despite of the absence 
of the time periodicity. However, we cannot now present any refine mathematical proof on 
quasibreather stability for infinite time. 

On the other hand, one often reveals that the lifetime of a quasibreather, if it actually 
decays in time, can be exponantionally large for small deviations from the exact breather. 
For example, in Fig. [TT], we depict the time evolution of the quasibreather in the weakly 
coupled Duffing chain which lives, at least, up to t ~ 10 6 . 

Finally, since we cannot tune exactly on the strictly time periodic breather solution in 
any physical experiment (and even in numerical experiments), quasibreathers seem to be 
much more relevant spatially localized dynamical objects than the exact discrete breathers. 
Some additional detail for this point of view can be found in 
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It is easy to prove numerically that, for resonable values of the parameters 7 and ol\ = x'i(O), 
Eqs. (|14M15p allow the breather solution, which is very close to that of the exact equations (|1Q|) . 
It is possible owing to the specific dependence of uj on the vibrational amplitude for the Duffing 
oscillator (see Fig. [6j) 
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